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ABSTRACT 

We present mid infrared (MIR) spectra of the Seyfert 2 (Sy 2) galaxy NGC 1808, 
obtained with the Gemini's Thermal-Region Camera Spectrograph (T-RcCS) at a spa- 
tial resolution of ~26pc. The high spatial resolution allowed us to detect bright poly- 
cyclic aromatic hydrocarbons (PAHs) emissions at 8.6^m and 11.3/im in the galaxy 
centre (^26 pc) up to a radius of 70 pc from the nucleus. The spectra also present 
[Ne n]12.8/im ionic lines, and H2 S (2)12.27^m molecular gas line. We found that the 
PAHs profiles arc similar to Peeters's A class, with the line peak shifted towards the 
blue. The differences in the PAH line profiles also suggests that the molecules in the 
region located 26 pc NE of the nucleus are more in the neutral than in the ionised state, 
while at 26 pc SW of the nucleus, the molecules are mainly in ionised state. After re- 
moval of the underlying galaxy contribution, the nuclear spectrum can be represented 
by a Nenkova's clumpy torus model, indicating that the nucleus of NGC 1808 hosts a 
dusty toroidal structure with an angular cloud distribution of a — 70°, observer's view 
angle i — 90°, and an outer radius of Ro ^0.55 pc. The derived column density along 
the line of sight is N# = 1.5 x 10 24 cm 2 , which is sufficient to block the hard radiation 
from the active nucleus, and would explain the presence of PAH molecules near to the 
NGC 1808's active nucleus. 
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1 INTRODUCTION 

Mid-infrared (MIR) spectra of active galactic nuclei (AGN) 
are dominated by the emission features at 3.3, 6.2, 7.7, 8.6, 
11.3 and 12.7/mi, c ommonly cal l ed th e unidentified infrared 
(UIR) bands (e.g. iGillett et all 1 19731 : iGeballe et all 1 19851 : 
ICohen et al.l 1198a ). and generally attribut ed to Polycyclic 
Aromatic Hydrocarbon (PAHs) molecules jLeger fc Puget I 
11984 iPuget fc Leger~lll989l : lAllamandola et al.ll 19891 . 1 19991 ). 
The spectra of AGN also present prom inent molecular hy 



drogen and forb i dden emission lines dSturm et al. 200d: 



Weedman et al. I 120051 : IWu et al. I 120091 : iGallimore et all 
2010; Sal es et al. |2010| ). as well as silicate bands at ~ 9.7 /im 
and 18 /im, both in emission and/or absorption. 

The absence of PAH emission close to the central engine 
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of active galaxies has been attributed to the destructive ef- 
fects of highly energetic UV photons (his ~ 40 eV) supplied 
by the active nucleus and, therefore, the presence of these 
features have be en mostly used as a t racer of massive star 
forma tion (e.g. iPeeters et al. I |2004|; ISiebenmorgen et al. I 



lAllamandola et al 
(|l992al lbl) has shown that 



1999; iTielens I 1200 



However, 
survive if the 



PAHs can 

hard AGN radiation is shielded by dusty g as with high col- 
umn density. Recently, ISales et al. I (|201Cf ) showed that a 
large fraction (> 70%) of Seyfert galaxies shows bright PAH 
emission lines, and demonstrated that PAHs in AGN present 
higher ionisation fraction and larger molecules (> 180 car- 
bon atoms) than in Starburst galaxies. 

The unified model of AGNs proposes the existence of a 
dense concentration of absorbing material surrounding the 
central engine in a toroidal distribution, which blocks the 
broad line region (BLR) from the line of sight in Seyfert 
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2 (Sy2) galaxies (see lAntonucci I 1 19931 s ) . Nowadays such 
toroidal struct ure is believed to be composed by gas clu mps 
or clouds (e.g. iHonig et al. 1120061; iNenkova et al. 1120021 ). In 
fact, very recently we (|SaIes et al. Il201ll ) have shown that 
the nuclear spectral energy distribution of NGC3281, a 
Sy 2 galaxy, can be well described by clumpy torus mod- 
els with about 14 dusty clouds in the torus equatorial ra- 
dius. Similar results are also obta i ned for other AGN (e.g. 
iRamos Almeida et al. Il2009l . 1201 ll ; iNikutta et aL~ll2009l ). 

It has been known by approximately five decades that 
NGC 1808 presents a peculiar nuclei, evidenced by the detec- 
tion of star forming regi ons disposed in a circumnuclear ring, 
the so called hot spots (ISersic fc Pastoriz al l 19651; I Pastoriza 
ll967l:J Sersic fc Pastoriza 1119671 ). In addition. lGalliano et al. 
1 20051 ) detected bright MIR sources in the central region 



HD26967 
NGC 1808 



of NGC 1808, assoc iated with already known radio sources 
(|Saikia et al. Ill99(il ). suggesting that these sources are young 
embedded star clusters. The presence of an active nucleus in 
this source is supported by the prese nce of extended wings 
on bo th H a and [N ll] emission lines (|Veron-Cettv fc Veronl 
1 19851 ). These authors found that NGC 18 08 falls in the 
Seyfert region of the iBaldwin et"al] (|l98lT) diagram with 
Ha/ [Nil] A 6583 = 0.94 pm 0.05 

and presents the full width at half maximum (FWHM) 
of Ha component approximately 550km/s. In addition, 
lAwaki et al. I (| 19961 ) have detected hard X-ray emission vari- 
ability suggesting the presence of a low-luminosity AGN in 
NGC 1808. This object is a nearby galaxy (d=13Mpc, us- 
ing Ho = 74km s" 1 Mpc -1 ) exhibiting evidence that strong 
star formation and an active nucleus coexist in the central 
regions, making it a key object to investigate PAH emission 
and torus structure through N band spectroscopy. 

Here we present a detailed study of N band spectra 
of NGC 1808 using ground-based MIR high angular reso- 
lution, R~ 100 spectra obtaine d with the Thermal- Region 
Camera Spectrograph (T-ReCS; lTelesco et al. Ill998l ) at the 
8.0m Gemini South telescope, achieving a spatial resolution 
of 0'.'36, or ~ 26pc at NGC1808 distance, which makes our 
data quite adequate to study the molecular and dust dis- 
tribution in the inner ~ 100 pc of this galaxy. This paper is 
structured as follows: in Section [2] we describe the observa- 
tion and data reduction processes. Results are discussed in 
Section[31 Summary and conclusion are presented in Section 

m 



2 OBSERVATIONS AND DATA REDUCTION 

The A^-band spectra of NGC 1808 were obtained with the 
T-ReCS in queue mode at Gemini South, in 2009 August 
28 and 29 UT, as part of program GS-2009B-Q-19. The ob- 
servational conditions were photometric, with water vapour 
column in the range 5-8mm, and image quality of 0'.'34 in 
the Af-band, measured from the telluric standard acquisition 
images observed before NGC 1808. The A-band luminosity 
profile of the galaxy is slightly extended with respect to that 
of HD 26967 telluric star (see Fig. [TJ 

We used a standard chop-nod technique to remove time- 
variable sky background, telescope thermal emission, and 
the effect of 1// noise from the array/electronics. In order 




Figure 1. The TV-band spatial emission profile (solid line) of 
NGC 1808 compared with that of the telluric star HD 26967 (dot- 
ted line). The fluxes were normalised to the peak value. 



to include only the signal of the guided beam position in the 
frame and avoid possible nod effects in the spatial direction 
the slit were oriented along P. A. = 45°with a chop throw of 
15", oriented along P.A. =45°. The same slit position/nod 
orientations was used for the telluric standards. The T-ReCS 
configuration used was the low-resolution (R ~ 100) grat- 
ing and the 0/36 slit, for a dispersion of 0.0223 pixel -1 
and a spectral resolution of 0.08 fim. The pixel scale is 
0.089"pixeP 1 in the spatial direction, and the slit is 2l'.'6 
long. The resulting spectral coverage is between 8- 13pm cen- 
tred at 10.5/im. The total on-source integration time was 
900s. 

We used the midir and GNIRS sub-packages of Gem- 
ini IRAF0 package to reduce the data. In addition, to ex- 
tract the final spectrum we combined the chop- and nod- 
subtracted frames using the tasks tprepare and mistack. 
Wavelength calibration was obtained from the skylines in 
the raw frames. To remove the telluric absorption lines from 
the galaxy spectrum, we divided the extracted spectrum for 
each obser ving night by that of the telluric standard star 
HD 26967 (|Cohen et al. Ill999l ). observed before the science 
target. The flux calibration used the task mstelluric that 
interpolates a blackbody function to the spectrum of the 
telluric standard in order to derive the proper instrumental 
response. 

In Fig. [2] we show the T-ReCS slit position super- 
imposed on our N band acquisition image, as well as 
the hard (> 2keV) and s o ft (< 1.5keV) X-ray images 
from iJimenez- Bailon et al~l (|2005l ). We can see that the 
X-ray images show two point-like emission sources off- 
set from the NGC 1808 nucleus (see in Figure |2j). These 
emitting regio ns were also detected in the radio and in - 
frared images (|Saikia et al. I ll990l ; iKotilainen et al. I Il996r i 
and have been attributed to supernova remnants (SI) 



1 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA), Inc., under cooperative agree- 
ment with the National Science Foundation. 
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and Hll regions (S2) (ISaikia et al. I fl9Q0l ; iKotilainen et al. I 
Il996l ; |jimenez-Bail6n et al. 120051 ). The T-ReCS slit position 
crosses the soft X-ray source SI and diffuse X-ray emission. 



gradient from the centre to the outer regions. The similarity 
between the luminosity profile of the PAH emission bands, 
ionic line and continuum suggests that the same mechanism 
dominates the excitation of both ionic and PAH bands. 



3 RESULTS AND DISCUSSION 

3.1 Emission Lines Measurements and Radial 
Gradients 



We extracted 7 one dimensional spectra spaced four pixels 
or 0//36 (~ 26 pc) along the slit: one centred in the unre- 
solved nucleus, and six extractions centred at 26 pc, 52 pc, 
and 78 pc to the northeast (NE) and southwest (SW) di- 
rections. The spectra are presented in Fig. [31 and clearly 
show intense emission in the PAH bands (8.6 fim, 11.3 pm 
and 12.7 /an), [Nell] 12.8 pm ionic line, as well as the molec- 
ular hydrogen rotational line, H2 S (2) at 12.27 /im, in the 
nuclear extraction, and at 26 pc in both NE and SW di- 
rections. Weak PAH bands are also observed at 52 pc SW 
direction of NGC 1808, while the outer extractions at 78 pc 
NE and SW do not show any emission lines. This result in- 
dicates that in this galaxy the PAH bands are emitted inside 
a region less than 70 pc from the active nucleus, and clearly 
show one of the most rich PAH features at 8.6, 11.3 and 
12.7 pm observed using T-ReCS so far in AGNs. 

To disentangle the different spectral components in each 
ID spectra, both emission lines and continuum, and deter- 
mine the contribution of each component to the spectral 
energy distribution (SED) of NGC 1808 as a function of the 
distance to the gal axy centre, we used the PAHFIT0 code 
jSmith et al. I [20071 ) routines. The input flux uncertainties 
required by PAHFIT were assumed to be 10% of the flux ob- 
serv ed, which are the expec t ed ones for T-ReCS observations 
(see iRadomski et al. 1120021 ; iMason et al~ll2006l ) . 

The pahfit code assumes that the spectrum is formed 
by continuum emission from dust and starlight, prominent 
emission lines, individual and blended PAH emission bands, 
and that the light is attenuated by extinctio n due to silicate 
;rains . The code use the dust opacity law of iKemper et al."l 
2004]), where the infrared extinction is considered as a 
power law plus silicate fe a tures peaking at 9.7 /im. For more 
details see ISmith et al. ] (|2007h . The r emaining input pa - 
rameters are the same as those used bv lSales et al. I (|2010t ). 
which are appropriate for AGNs. The results of the decom- 
position of the NGC 1808 spectra are shown in Fig. f4j and 
the derived emission line fluxes are given in Tabs.[T]and[2] It 
is worth note that 10.7pm PAH emission at the centre and 
26 pc NE is at the detection limits and its emission profile 
can not be seen in the PAHFIT fitting at Fig. [4] In addition, 
note that the 11.3pm PAH emission complex is composed 
by 11.2pm and 11.3pm emission bands, which can be seen 
in Fig. El 

Fig-Elshows the PAH bands and [Ne 11] 12.8 pm emission 
line integrated fluxes, as well as the 10pm flux continuum 
emission flux, as function of distance to the nucleus. Both 
emission lines and continuum present a strong decreasing 



3.2 NGC 1808 PAH Profiles Analysis 

iPeeters et all (|2002T l and Ivan Diedenhoven et al~l (|2004l ) 

demonstrated that the PAH bands in a sample of reflec- 
tion nebulae, Hll regions, young star objects, evolved stars 
and galaxies can show different profiles, which they pro- 
ceeded to classify based on their shape and peak posi- 
tion into three classes. The A class has the peak shifted 
to the blue, while B and C would show the peak cen- 
tred and shifted to the red, respectively. According to these 
authors, Hn regions, non-isolated Herbig AeBe stars and 
young stellar objects generally can be classified as A class 
for all 3 — 12pm PAH emission features. Their galaxy sam- 
ple shows the same profile as the Hn regions (A class) ex- 
cept for the 11.3pm profile, which resembles that of some 
evolved stars (B class). They also point out that most 
planetary nebulae and post-asymptotic giant branch stars 
belong to the second group (B class) and just two pecu- 
liar pos t-asymptotic giant branch stars form a third clas s 
(C class IPeeters et ah 20021; van Diedenhoven et al. 2004 ) . 
Howev er JPeeters et al. I (|2002T l and lvan Diedenhoven et al. I 
|2004l ) do not analyse the PAH line profiles in the AGNs 
environment. These authors have been attributed the PAH 
different profiles to differences in the physical conditions, 
composition of the PAH family and/or the accumulated ef- 
fect of process ing in the regions where the emission origi 
nates (see alsolGenzel et al. Ill998l;lvan Diedenhoven et al. 



| 2004 l; iLebouteiller et al. 1 120071 ; iTielens I l20Qg| ; ISales et al 
[2010). In addition, other authors propose that the trends in 
the PAH band line profiles and intensities are a function of 
the PAH molecules size, charge, geometry, an d hete rogeneity 
|Honv et al. Il200ll ; IBauschlicher et al. ||200S| , 120091 ) . 

In order to investigate the PAH profiles using high spa- 
tial resolution (~26pc) of the Seyfert NGC 1808, we isolate 
the PAH emission subtracting a general and local continuum 
under th e 8.6 pm and 1 1 .3 nm band s, using the same proce- 
dures of IPeeters et al. I |2002l ) and Ivan Diedenhoven et al. I 
|2004l ). Figures [6] and [3 show the 8.6 pm and 11.3 pm PAH 
profiles from our dat a, compared to repres entative PAH pro- 
files of the A and B IPeeters et al~l (|2002h classed- 

In general, the 8.6pm and 11.3pm PAH emission line 
profiles in NGC 1808 resemble Peeters' A class. However, 
the 26 pc SW and 26 pc NE 8.6pm PAH profiles are wider 
than the representative A class example, and also shifted to 
shorter wavelengths. Compa r ison w ith the theoretical mod- 
els from IBauschlicher et al. 1 (|2009T l would suggest that the 
excess emission in the blue wing of the 8.6pm line profile is 
evidence that the off-nuclear regions at 26 pc SW and 26 pc 
NE contain more ionised PAH material than the nuclear re- 
gion in this galaxy. The same conclusion can be derived from 
the equivalent width of the 8.6pm assuming that this line is 



2 Source code and documentation for PAHFIT are available at 3 The IPeeters et al. ] ||2002|) profiles are available at 

http: / /tir. astro. utolcdo.edu/jdsmith/rcscarch/pahfit.php |http: / /www. astro. uwo.ca/~ep ectcrs/rcscarch-data,shtml 
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Figure 2. The top panels shows the T-ReCS long-slit position over-plotted on the contoured NGC 1808 acquisition image. The contours 
are linear and stepped at 20% of the peak. The bottom pan els show the T-ReCS slit superposed to the Chandra ACIS soft (left) and 
hard (right) X-ray images from [jimenez-Bailon et al. I (120051 ). 



Table 1. PAH Emission Line Strengths (10 13 Wm 2 ) and Equivalent Widths (/tm). 



Position 


Label 


8.3//m 




8.6/im 




10.7/1 


m 


11.3/1 


m 1 


12.0/i 


m 


12.7/1 


m 


Flux 


EW 


Flux 


EW 


Flux 


EW 


Flux 


EW 


Flux 


EW 


Flux 


EW 


26 pc SW 


-1 


3.25±0.1 


0.562 


3.15±0.1 


0.517 






2.95±0.2 


0.396 


1.19 ±0.1 


0.163 


2.05 ± 0.3 


0.291 


Center 


(] 


5.35±0.S 


0.119 


11.2±0.7 


0.242 


0.29±0.3 


0.006 


11.8 ± 1.2 


0.251 


4.85 ±0.7 


0.110 


7.54 ± 1.6 


0.182 


26 pc NE 


1 


4.33 ±0.3 


0.328 


4.44 ±0.2 


0.334 


0.11 ±0.1 


0.009 


5.70±0.4 


0.477 


2.46 ±0.5 


0.208 


3.46 ±1.1 


0.285 


52 pc NE 


2 






3.13±0.1 


0.990 






1.32±0.1 


0.775 


0.81 ±0.1 


0.493 


0.39±0.1 


0.265 



"'"The ll.3ji.rn PAH emission complex is composed by 11.2pm and 11.3/xm emission bands. 

Table 2. Ionic and Molecular Emission Line Strengths (1Q~~ 14 W m~ 2 ) and Equivalent Widths (/mi). 



Position 


Label 


[Arlll]8.99/jm 


[Siv]10 


5/im 


H 2 S(2)12.27/im 


[Ncl 


12.8/jm 


Flux 


EW 


Flux 


EW 


Flux 


EW 


Flux 


EW 


26 pc SW 


-1 


0.20 ±0.10 


0.004 


0.28±0.19 


0.005 


0.79±0.25 


0.016 


8.08 ±0.46 


0.170 


Center 


(] 






0.64±0.10 


0.001 


3.13 ±1.40 


0.011 


45.6±2.19 


0.166 


26 pc NE 


1 






0.28±0.11 


0.001 


0.98 ±0.46 


0.012 


17.5±9.17 


0.210 


52 pc NE 


2 


3.88±0.18 


0.235 


2.80±0.13 


0.238 


0.24 ±0.08 


0.022 


3.18±0.16 


0.326 
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Figure 3. Spectra of NGC 1808 extracted in 26 pc steps along P.A. = 45°. The dashed lines indicate the positions of the PAH bands, 
ionic as well as hydrogen molecular lines that arc listed at Tabs, [l] and [2] The telluric O3 band is represented by the hatched area. 



emitt ed by ionized PAH molecules (see lAUamandoIa et all 
Il999l ). 

Fig. [7]sliows that the 11.3/tm PAH profile from the 26 pc 
SW spectrum is broader towards short er wavelengths than 
the represen tative A class. According to lBauschlicher et al. I 
(|200Sl . l2009l ), the emission creating the shorter wavelengths 
of the 11.3/xm profile would be produced mainly by neutral 
PAH molecules, and we would conclude that the amount of 
neutral emitting PAH material in the 26 pc SW region is 
larger than in the NGC 1808 centre. The same conclusion 
can be derived from the equivalent width of the 11.3/im 
assu ming that this line is emi tted by neutral PAH molecules 
(see lAUamandoIa et alii 19991 ). 

In addition, we can see from the values in the Ta- 
bled] that the equivalent width (EW) of the PAH at 8 6/xm 
which is due to the ionized PAHs fe.g.|p~uget fc Leger II 19891; 
Allamandola et al.|[l989l . Il999l ; IPeeters et al. ll200Sj ; iTielens I 



2008? ) . is larger at 26 pc SW extraction, while the largest 



EW of the neutral PAHs (see references above)is at 26 pc 
NE extraction. In summary, the comparison of the observed 
PAH profiles and their EWs with the existing observational 



and theoretical studies indicate that at 26 pc SW of the 
nucleus in NGC 1808, the ratio of ionised to neutral PAH 
molecules is larger (1.3) than in the nuclear region (0.96). 
On the other hand, the EW ratio of ionised to neutral PAH is 
smallest, EW11.3/nm/EW8.6/Lim~0.7, at 26 pc NE suggest- 
ing that this region has more PAHs in the neutral ionization 
stage. Interestingly, the region hosti ng more ionised PAH 
coinci des with the ST X-ray source of lJimenez-Bailon et al. I 
i|2005f ). while the side hosting more neutral molecules coin- 
cides with the extended diffuse X-ray emission (see Fig. [2}. 
Overall, the asymmetry in the PAH profiles correlates with 
regions associated with soft X-ray emission (E < 1.5keV). 

Our data make it quite evident that the 8.6/im and 
11.3/mi PAH emission (from ionised and neutral molecules) 
are present in regions of soft and hard X-ray emission, as well 
as near to the active nucleus itself (r <30pc). This immedi- 
ately implies that both ionised and neutral PAH molecules 
are surviving in the hard X-ray radiation field of the AGN, 
as predicted if a dusty structure, such as the putative torus 
in the AGN unified model, is present to partially block the 
strong radiation field from the central black hole. In this 
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8 9 10 11 12 8 9 10 11 12 13 

Wavelength (jum) Wavelength (fim) 

(c) (d) 

Figure 4. Result of the PAHFIT decomposition of the NGC 1808 spectra: (a) 26 pc SW; (b) unresolved nucleus; (c) 26 pc NE: and (d) and 
52 pc NE extractions, respectively. The data points are represented by the dots, with uncertainties plotted as vertical error-bars. The 
best fit model is represented by the blue line. The dotted black line indicates the mixed extinction components, while solid green and 
orange lines represent total and individual thermal dust continuum components, respectively. The violet and magenta lines represent the 
ionic and PAH emission lines. Note that PAHFIT code uses the violet solid line to represent the H2 emission as well. 



context, it is an interesting issue to investigate if the nu- 
clear spectrum of NGC 1808 present evidence of dusty torus 
that would explain the presence of PAH molecules near to 
the active nucleus. 



3.3 Is there a Dusty Torus Present in the Central 
Region of NGC 1808? 

The presence of a nuclear dusty toroidal structure, mainly 
composed by silicate and graphite grains, leaves un- 
mistakable signatures in the observed SEDs of AGN. 
While the sublimation of the graphite grains creates IR 
emission at A ^ 1/zm, the ~9.7/im feature observed in 



Barvainis I 119871: IPier fe Krolik \ 1 19921 ; 
1994 ISiebenmorgen et al. I I2005T 



Rodriguez- Ardila fc Mazzalavl 120061 ; 



Granato & Danese 1 


Fritz et al. 




120061; 


Riffel et al 




20091). 



Some authors d escribe the torus as a continuous density 
distribution fe.g.lPier fc Krolik Hl992l;lGranato et al. Ill997l ; 



been postulated that, for dust grains to survive in the torus 
environment, they should be formed in clumpy structures 
|Krolik fc Begelmanl Il988h and once pro vide a natural 
attenuation o f the silicate feat ure (e.g. iNenkova et al. I 
12001 l2008al lbl iHonig et al 1120061) . 

Sales et al. I (|201ll ). we removed the under- 



Following 



lying host galaxy contribution from the NGC 1808 nuclear 
spectrum by subtracting the average spectrum of the two ad- 
jacent extractions (apertures -1 and 1 in Fig. [5J. Moreover, 
we masked the emission lines and the telluric band region 
(Fig[3J) using a simple interpolation, and then compared the 
resultin g pure nuclear spectrum agai nst the Nenkova's mod- 
els (e.g. INenkova et al. |[20o3 . l2008al fbh . 

These modelsQ assume that the torus is formed by dusty 
clumps constrained by the following parameters: (i) the 
number of clouds, No, in the torus equatorial radius; (ii) the 
optical depth of each cloud in the V band, tv\ (iii) the radial 



distribution fe.g.lPier fc Krolik 111992k iGranato et al. II 19971 ; 4 The models are 

ISiebenmorgen et al. 2004 : Fritz et al. 20061 ) . but it has |http: / /www.pa.uky.edu/clumpy/| 



available 



from 
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Figure 5. Emission lines and lO^tm continuum luminosity profiles along the slit. 
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Rest Wavelength (um) 

Figure 6. Comparison between the NGC 1808 8.6/im PAH profiles and the As. 6, and Bg.g Peeters's PAH classes. All profiles have been 
normalized to the peak intensity. The solid black line correspond to the observed PAH profile. 



extension of the clumpy distribution, Y — Ro/Rd, where Ro 
and Rd are the outer and inner radius of the torus, respec- 
tively; (iv) the radial distribution of clouds as described by 
a power law oc r~ q ; (v) the torus angular width, constrained 
by a Gaussian angular distribution described by a half-width 
<r; («) the observer's viewing angle i. 

The ionic and PAH emission lines, as well as the telluric 
band (Fig[3]) were ignored for the fit. The final spectrum was 
compared with the ~ 10 6 clumpy theoretical SEDs, and 
the best fit is obtained by searching for the minimum in the 
equation: 



-T 



N ' W -I? ^ 2 

-Fobs, Aj ^mod, A; 



1=1 v z 



(1) 



where N is the number of data points in the spectrum, 



Fobs, \i and F mo£ j, x 4 are the observed and theoretical 
fluxes at each wavelength, and S\ i are the uncertainties 
in Fobs, Xi , calculated as 10% of F\ (see Radomski et al.1 
120021 ; iMason et al~l 120061 ; ISales et al. I I2OIII ). Both F nh , >~ 
and Fmod.Xf were normalised to unit at 11.0/im, with the 
uncertainties correctly propagated. Table[3]show the param- 
eters of the best fit, and Fig.|5]shows the resulting theoretical 
SED superimposed on the NGC 1808 pure nuclear spectrum. 

The mean and standard deviation of the model param- 
eters were calculated from the subset of models with a x 2 
value within 10% of the best-fitting result (144 in total). 
The locus of the corresponding theoretical SEDs is plotted 
as a gr ey region in Figure [8] We used a similar approach to 
that of lNikutta et al. I (|2009l ), who test 5%, 10%, 15% and 
20% x 2 deviation fractions, and found that our methodol- 
ogy results in a narrower locus of acceptable solutions. The 
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Figure 7. Comparison between the NGC 1808 11.3/im PAH profiles and the An. 3, and Bn.3 Peeters's PAH classes. All profiles have 
been normalized to the peak intensity. 




ISales et al. I 1201 ll ). This value is a factor of ~6 larger 
than that inferre d from X-ray observation (Lx-ray = 
1.6 x 10 40 erg s _1 Ijimenez-Bailon et al. 20051) afte r appl y 
the correction factor of 20 sug gested by I Elvis* et al. 1 1|2004) . 
iRamos Almeida et al. I (|2009l ) have also found similar re- 
sults from their clumpy torus model fitting (iJx-ray = 
6.8 x 10 41 erg s" 1 ) of the NGC1808 photometric data. How- 
ever, such discrepancy could be at tributed to the atte nua- 
tion from circumnuclear dust (e.g. iKrabbe et all 119941 s ) in- 
cluded in the la rger (850 pc) aperture used for the X-ray 
(850 pc) used by[j lmcnez- Bailon et al~1 (|2005l V 

Note that the high column density derived from the 
models would be e nough to blo ck the hard radiation from 
the central source (|Voit 1992a b), explaining how the ob- 
served PAH molecules can survive so close (r<26pc) from 
the Seyfert nucleus. 



Figure 8. Comparison between the NGC 1808 nuclear spectrum 
after removal of the underlying galaxy contribution (solid black 
line), and the best fit torus model (red line) as defined in the text, 
normalised at 11.0/im. The grey dotted lines correspond to the 
10% of best-fitting SEDs. 



mean values for the model parameters within the adopted 
X 2 deviations are shown in Table [3] 

The parameters obtained for the best fit model suggest 
that NGC 1808 hosts a dusty toroidal structure, where the 
dusty clouds, with individual optical depth of ry = 100 mag, 
occupy a toroidal volume of Ro/Rd = 10. The density pro- 
file distribution of the clouds follows a power-law r~ 1,5 , and 
there are 15 clouds in the equatorial radius of this torus. 
The toroidal structure is constrained in angular width by a 
Gaussian distribution of width a = 70°. 

According to our best-fit model, we would be look- 
ing along the torus equator (i = 90°), resulting in a 
line of sight column density of Nh ~ 1-5 x 10 24 cm~ 2 
and torus outer radius is Rn ~0.55 pc. We can also esti- 
mate the bolometric luminosity fro m clumpy torus model 



as Li, 



,9xl0 42 erg s" 1 (see iNenkova et al. I [2008bl ; 



4 SUMMARY AND CONCLUSIONS 

In this paper we analyse N band spectra of the Seyfert 2 
NGC 1808, obtained with the Gemini T-ReCS at a spatial 
resolution ~26pc. The central regions of this galaxy show 
extended PAH band emission, and conspicuous differences 
in the PAH line profiles are detected in scales of dozens of 
parsecs. Our main conclusions are: 

(i) The spectra of the unresolved nucleus and of re- 
gions located 26 pc in both NE and SW directions along 
the slit show intense 8.6/im, 11.3/im and 12.7/im PAH 
emission bands. The spectra also show [Ne 11] 12.8/zm and 
H2 S (2)12.27^iin emission lines. Both ionic and molecular 
emission are extended up to 70 pc from the nucleus. 

(ii) The PAH, [Ne n] 12.8/^m, and 10/im continuum emis- 
sion decrease radially in intensity. The similarity in their lu- 
minosity profiles would indicate that their excitation mech- 
anism is the same. 

(iii) The analysis of the PAH emission lines shows that 
the 8.6^im and 11.3jum line profiles in the central region 
are similar to Peeters's A class. However, the profile of the 
lines originating from the 26 pc SW and 26 pc NE regions 
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Table 3. Resulting Parameters for the fit of the CLUMPY models to the nuclear spectrum of NGC 1808. 



Parameter 




Best Pit 


Average 


Average 


Average 


Average 








5% 


10% 


15% 


20% 


Torus angular width (cr) 




70° 


65° ±5° 


62° ± 7° 


61° ±8° 


61° ±8° 


Radial extent of the torus (Y") 




1(1 


10 ± 


11 ± 7 


11±8 


11 ± 10 


Number of clouds in the torus equatorial rad 


ius (No) 


16 


14 ± 1 


13 ± 1 


13 ± 2 


12 ±2 


Power-law index of the radial density profile 


(?) 


1.5 


1.5 ± 0.1 


1.1 ±0.3 


0.7±0.5 


0.6±0.5 


Observers' viewing angle (i) 




90° 


80° ±9° 


76° ± 11° 


74° ± 13° 


71° ± 15° 


Optical depth for individual cloud ("Ty) 




100 mag 


100 ± mag 


99 ± 6 mag 


99 ± 6 mag 


100 ± 7 mag 


Total number of solution 






25 


144 


343 


525 



are different, implying that in the ratio of ionised to neutral 
PAH molecules in the 26 pc NE region is higher than in the 
26 pc SW one. 

(iv) Strong PAH molecular emission is present in the nu- 
clear extraction, as well as in the off-nuclear regions coinci- 
dent with hard and soft X-ray emission. This suggests that 
the PAH molecules in the nucleus of NGC 1808 are able to 
survive the AGN hard radiation field. 

(v) The nuclear N band spectrum after subtraction of 
the underlying galaxy contribution can be r epresented by 
a clumpy torus model (|Nenkova et al. Il2008bl ) . The torus in 
NGC 1808 would present a toroidal volume of Ro/Rd = 10, 
with each cloud having an optical depth tv = 100 mag and 
the distribution of the clouds following a r -1 ' 5 power-law. 

(vi) According to the best-fit torus model, our line of sight 
would be along the torus equatorial radius (i — 90°), result- 
ing a hydrogen column density of Nh ~ 1.5 x 10 24 cm" 2 . 
The torus would have an outer radius of Ro ~0.55 pc. Our 
modelling is consistent with NGC 1808 hosting a Sy 2 nu- 
cleus. 
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